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The  dynamics  of  electrons  in  a  low-frequency  wave  propagating 
perpendicular  to  a  uniform  magnetic  field  are  studied  and  the  impli¬ 
cations  of  these  results  for  transport  and  heating  by  the  lower 
hybrid  drift  instability  are  explored.  Below  a  threshold  ^ /T^  ^  y 
.25-. 5,  all  electron  energy  and  momentum  exchange  with  the  wave  are 
reversible  and  no  plasma  transport  is  possible.  Above  this 
threshold,  trapping  of  electrons  by  the  wave  potentials  takes  place 
and  causes  irreversible  electron  heating  and  momentum  exchange. 

These  results  imply  that  anomolous  transport  in  inhomogeneous 
plasma  with  weak  drifts  (diamagnetic  velocity  less  than  the  ion 
thermal  velocity)  may  be  substantially  less  than  previously  predicted. 
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1 .  INTRODUCTION 


Jhe  subject  of  heating  and  transport  in  high-temperature,  magnetically- 

confined  plasma,  where  classical  collisional  processes  are  weak,  is  of  great 

importance  both  in  achieving  thermonuclear  fusion  in  laboratory  confinement 

systems  and  in  understanding  fundamental  natural  phenomena  such  as  magnetic 

substorms.  Instabilities  whose  source  of  free  energy  is  the  plasma  density 

or  magnetic  field  gradients  are  likely  to  produce  fluctuating  electric  or 

magnetic  fields  which  cause  crossfield  particle  transport.  In  systems 

1/2 

with  rather  sharp  density  gradients  where  p./L  >  (m  /m.)  , 

-  I  n  e  1 

and  being  the  ion  Larmor  radius  and  density  scale  length,  respectively, 

the  lower-hybrid-drift  instability  is  unstable  and  is  expected  to  play  a 

dominant  role  in  the  evolution  of  the  plasma  profiles.^  "*  Such  sharp  gradients 

occur  in  a  variety  of  physical  systems  including  laboratory  plasmas  such  as 

9-pinches  and  reversed-f ield  pinches  and  space  plasmas  in  the  bow  and  tail 

of  the  earth's  magnetosphere.  The  linear  properties  of  this  mode  have  been 

1-3 

extensively  investigated.  The  instability  is  driven  by  the  pressure  gradient 

and  is  characterized  bv  a  frequency  w  and  growth  rate  y  given  by  u  %  kv^ 

%  >:>  Y  <  w  ,  where  Vj..  is  the  ion  diamagnetic  velocity,  is 

the  ion  gyrof requency  and  is  the  lower  hybrid  frequency.  The  growth  rate 

of  the  instability  is  sharply  peaked  perpendicular  to  ^  (k(|  =  £  0)  since 

electron  Landau  damping  for  finite  k |(  is  strongly  stabilizing. 

The  transport  and  heating  associated  with  the  nonlinear  evolution  of  the 

lower  hvbrid  drift  instability,  which  can  be  described  bv  an  effective  "anomolous" 
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resistivity,  have  been  calculated  in  a  quasilinear  analysis.  Particle  simu¬ 
lations  of  the  instability  in  both  straight  (9-pinch)  and  reversed  magnetic 
fields  have  also  demonstrated  that  electron  heating  and  transport  can  ac¬ 
company  the  nonlinear  development  of  the  mode.^’^ 


Heating  and  transport  are,  by  definition,  irreversible  processes  and  the 
subject  of  irreversibility  must  be  addressed  in  discussing  any  transport  or 
heating  mechanism.  In  low  temperature  plasmc  collisions  guarantee  irreversi¬ 
bility  on  the  macroscopic  scale.  In  collisionless  plasma,  resonant  wave 
particle  interactions  lead  to  irreversible  exchanges  of  energy  or  momentum 
(and  consequently  transport).  In  the  previous  quasilinear  investigation  of 

heating  by  the  lower  hybrid  drift  mode,  only  perpendicularly  propagating  waves 

2 

(k|(  =  0)  were  considered.  Ions,  which  behave  as  if  they  are  completely  un¬ 
magnetized  (since  w  >>  £2^),  resonantly  interact  with  the  wave  and  exchange 
momentum  and  energy.  The  electrons,  however,  which  are  tightly  bound  to  the 
magnetic  field  lines  (w  <<  £2g)  are  nonresonant  and  therefore  undergo  no 
irreversible  energy  or  momentum  exchange  (neglecting  7B  resonances).  The 
quasilinear  electron  "heating"  previously  calculated  simply  results  from  the 
coherent  sloshing  of  the  electron  distribution  function  in  the  lower  hybrid 
waves  and  is  completely  reversible.  In  particular,  for  nonresonant  particles, 
the  heating  rate  is  proportional  to  9 | E |  /3t,  and  the  temperature  increase 

■  'Xi  |  2 

scales  as  |E|  ,  E  being  the  amplitude  of  the  electric  field  perturbation. 

After  an  entire  growth  and  damping  cycle  of  the  wave  is  completed,  |e|  -+■  0 

and  there  is  no  net  electron  heating. 

Although  the  quasilinear  theory  predicts  that  no  irreversible  electron 

heating  can  take  place  for  Jc  •  ^  =  0,  computer  simulations  of  the  lower  hybrid- 
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drift-instability  indicate  that  electron  heating  does  actually  occur. 

These  simulations  are  carried  out  in  the  two-space  dimensions  orthogonal  to  B 
and  consequentlv  k|(  =  0.  If  this  heating  is  a  real  physical  effect  and  not  a 
consequence  of  numerical  errors.  It  must  be  associated  with  strong  non linearities 
in  the  electron  motion  not  property  described  by  the  quasilinear  theory. 

To  develop  an  understanding  of  electron  heating  by  the  1 ower-hvbr id-drif t 
instability,  we  study  the  electron  dvnnmics  in  a  single  i arge-amplitude,  low 
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frequency  (ui  <<  ftg)  wave  propagating  perpendicular  to 


Ey(y»t)  =  Ey(t)  cos(kyy  -  wt)  , 


where  the  time  dependence  of  Ey(t)  represents  the  growth  or  damping  of  the  wave. 
The  single  mode  approximation  is  a  reasonable  model  of  the  simulations  which 
were  largely  dominated  by  a  single,  coherent  wave.  Other  physical  effects 
including  VB  drifts,  the  two  dimensionality  of  wave  and  the  self-consistent 
evolution  of  the  wave  amplitude  are  neglected. 

At  small  wave  amplitude,  the  electron  motion  in  this  electric  field  is 
simply  given  by  the  usual  E  x  B  and  polarization  drifts  and  is  accurately 
described  by  the  quasilinear  theory.  Above  a  threshold  given  approximately 


k  (cE  /Bn  )  =  k  Av  >  1  , 

y  y  e  y 


the  electron  cyclotron  motion  is  strongly  modified  and  the  electron  motion 
becomes  stochastic.  The  distance  Ay  in  Eq.  (2)  is  the  y  displacement  of  the 
electron  due  to  the  polarization  drift.  When  this  displacement  is  comparable 
to  the  wavelength  of  the  electric  field,  the  strong  modification  of  the 
electron  dynamics  should  not  be  too  surprising.  Above  the  threshold  given  in 
Eq.  (2)  substantial  electron  heating  can  occur. 

The  present  investigation  of  electron  dynamics  is  related  to  previous 

studies  of  stochastic  ion  motion  in  large  amplitude  waves.  Smith  and  Kaufman 

studied  the  ion  transition  to  stochastic  motion  in  an  obliquely  propagating 
8 

wave.  Stochastic  particle  motion  in  perpendicularly  propagating  waves  has 
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also  been  investigated.  However,  the  authors  were  primarily  interested  in 
lower  hybrid  heating  of  ions  and  therefore,  only  considered  cases  where  w  >> 

More  recently,  heating  by  a  large  amplitude  standing  wave  has  been  studied. *^* 

The  present  calculation  is  actually  somewhat  simpler  than  the  previous  in¬ 
vestigations.  The  existence  of  two  disparate  time  scales  ^e*  and  w  *  allows 
a  straightforward  analytic  investigation  of  the  electron  motion  and  subsequently 
to  the  threshold  given  in  Eq.  (2). 

In  Section  II  of  this  paper,  the  electron  motion  and  transition  to  stoch- 
asticity  in  a  single  wave  is  discussed.  In  Section  III,  these  results  are 
extended  to  an  arbitrary  number  of  waves  in  one-space  dimension.  In  contrast 
with  some  turbulence  theories  based  on  the  renormalization  of  particle  orbits**, 
spatial  diffusion  is,  strictly  speaking,  not  possible  in  one-dimensional,  low- 
frequency  turbulence.  In  Section  IV,  the  implications  of  these  results  for 
the  saturation  and  transport  associated  with  the  lower-hybrid-drift  unstability 
are  discussed.  In  Section  V,  the  essential  results  and  conclusions  are  summarized. 
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II.  STOCHASTIC  ELECTRON  MOTION  IN  A  SINGLE  WAVE 


The  equations  of  motion  for  electrons  in  a  plane  wave  propagating 
pendicular  to  a  uniform  magnetic  field  B  =  Bz  are  simply 

y  =  ft  x  -  eE  /m  (3) 

J  e  y  e 

X  =  -S!  v  (4) 

e- 

where  ^e=eB/mec  is  the  electron  gvro frequency,  E^Cv.t)  is  given  in  Eq.  (1), 
and  the  dot  denotes  d/dt.  Equation  (4)  can  be  integrated  once  and  x  can  then 
be  eliminated  from  Eq.  (3)  to  obtain  the  single  equation 

0  +  9  *  a  cos(0  -  vt  +  <J>)  (5) 

where  0  =  k  y,  v=a)/s l  ,  a  =  -k  eE°/mf!^,  time  has  been  normalized  to  the 

y-'  g  y  y  0 

gyrofrequency  and  <t>  is  the  initial  phase.  This  equation  has  been  studied 

9 

previously  in  the  limit  v  >>  1.  We  consider  the  opposite  limit  v  <<  1. 

In  the  limit  v  -+  0,  Eq.  (5)  can  be  integrated  exactly  to  obtain  the 
particle  energy 


r)  n 

H=0  /2+9  /2-a  sin(0  +  $)  =  const.  .  (6) 

The  particle  motion  can  be  simply  understood  by  plotting  the  constant  H 
curves  in  the  0-0  phase  plane  as  shown  in  Fig.  la  for  a  <<  1  and  Fig.  lb  for 
a  >>  1.  The  electrons  simply  move  along  the  constant  H  curves  in  this  phase 
space.  When  a  <<  1,  these  curves  are  essentially  concentric  circles  weakly 
modified  by  the  presence  of  the  wave.  The  electric  motion  around  these  curves 
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on  a  time  t  1  corresponds  to  the  usual  electron  Larmor  motion.  When 
a  >>  1,  the  phase  space  structure  becomes  more  complex  with  the  formation 
of  x  and  0  points  as  can  be  seen  in  Fie.  lb.  The  electrons  aeain  circulate 
around  the  closed  phase  space  curves  on  a  time  t  ^  1 ,  except  in  narrow  bands 
around  the  separations  where  the  period  becomes  very  long.  In  this  limit  the 
usual  Larmor  orbit  is  strongly  modified  by  the  presence  of  the  wave.  The 
islands  shown  in  Fig.  lb  correspond  to  regions  where  the  perpendicular 
electric  field  is  large  enough  to  overcome  the  magnetic  field  and  trap  the 
electrons. 

The  location  of  the  stationary  points  (x  and  0  points)  in  the  phase  space 
of  Fig.  1  can  be  calculated  by  solving  3H/30  =  3H/30  =  0  or 

es  -  0  (7a) 

0  =  a  cos(0  +  <f>)  .  (7b) 

^  s 

The  solutions  to  Eq.  (7b)  are  shown  graphically  in  Fig.  2.  For  a  <  1,  there 
is  a  single  0  point  as  shown  in  Fig.  la  given  by 

0g  =  a  cosij)  <<  1  ,  (8) 

while  for  a  >  1,  there  are  multiple  x  and  0  points  increasing  in  number  with 
a.  For  a  >>  1,  the  solutions  are  approximately  given  by 

0  g  =  —  1 4>  +  (n  +  |)ti][1  +  (— l)n/ot]  ,  (9) 

where  n  is  an  integer  and  we  have  assumed  |0g|  a<<l.  The  outermost  0  point 
is  given  approximately  by  0g^ct  (see  Fig.  2). 


8 


The  electron  motion  in  the  vicinity  of  the  0-points  of  Fig.  lb  can  be 
investigated  by  expanding  H  in  a  Taylor  series. 


H  i  02/2  +  uil(0-6  )2/2  , 
b  s 


(10) 


222  22  1/2  1/2 

where  =  3  H/30  =  [l+aCl-S^/a  )  ]  is  the  bounce  frequency  of  the 

electron  in  the  potential  well.  For  1 8  |  «a  (large  islands  in  Fig.  lb), 

1/2  1/2 

oj^  -  (1+a)  =a  ,  which  corresponds  to  the  usual  bounce  frequency 
o  1/2 

ky(e<{i  /m^)  '  of  a  trapped  electron  in  the  absence  of  the  magnetic  field. 

Note  that  since  col,  the  bounce  frequency  for  these  electrons  is  greater  than 
the  gyrofrequency .  The  bounce  frequency  of  electrons  in  the  smaller  islands 
of  Fig.  lb  decreases  monotonicallv  with  increasing  | © ^ |  until  %  1  for  the 
outermost  islands. 

When  v^O,  the  particle  energy 


H  =  0 2 / 2  +  62/2  -  a  sin(9  -  vt  +  <t)  ,  (11) 


the  generalization  of  Eq.  (6)  for  finite  v,  is  no  longer  a  constant  of  the 
motion.  However,  in-so-far  as  v  <<  1 ,  the  energy  is  approximately  conserved 
during  the  rapid  gyro-motion  discussed  in  Fig.  1  and  this  rapid  motion  is 
basically  unchanged  by  v,  which  simply  enters  Eq.  (H)  as  a  phase  shift. 

On  a  longer  time  scale  vt  'v  1,  the  particle  energy  change  cannot  be  neglected. 
To  simplify  the  discussion  of  the  electron  dynamics  for  v  /  0,  the  two  cases 
o  >>  1  are  considered  separately. 

A.  Small  Amplitude  Waves:  q  <<  1 

With  v  ^  0  the  constant  H  curves  shown  in  Fig.  la  change  on  a  time  scale 
v  1  and,  in  particular,  the  0-point  oscillates  periodically 
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0g  ~  a  cos(4>  -  ut).  (12) 

The  particle  energy  H,  which  is  no  longer  constant  can  be  calculated  directly, 

H  =  -a  sin(0  -  vt  +  41)  +  va  cos(0  -  vt  +  <j>)  ,  (11) 

where  terms  not  proportional  to  a  or  v  have  been  eliminated  by  invoking 
Eq.  (5).  To  find  the  long  time  evolution  of  H,  we  average  Eq.  (13)  over  the 
fast  cyclotron  motion  contained  in  0(t), 

<H>t  =  -a  <sin(0  -  vt  +  a)>t  +  va<cos(0  -  vt  +  <j>)>t  .  (14) 

For  electrons  near  the  O-point  of  Fig.  la,  i.e.,  electrons  with  small  Larmor 
radius,  <sin(0  -  vt  +  <j>)>t  -  sin(eg  -  vt  +  < t>)  and 

<H>t  -  -a  sin(0g  -  vt  +  <(>)+  va  cos(0g  -  vt  +  <J>)  (13) 

=  d  H(0s,t)/dt  , 

where  H(9s,t)  is  the  energy  at  the  O-point.  Since  the  rate  of  change  of  the 
energy  of  the  particle  near  the  O-point  is  the  same  as  the  energy  of  the  0- 
point,  the  particle  must  follow  the  O-point.  As  the  O-point  in  Fig.  la 
oscillates  in  0,  the  particles  simply  follow  along. 

The  average  energy  gain  of  the  electrons  over  many  oscillations  v  ^  can 
be  obtained  by  further  averaging  Eq.  (15)  over  the  time  scale  v  ^  under  the 
assumption  that  a/a  <<  v.  We  find 

<<H>>t  -  -  aa/2  =  -  d  a^/dt  ,  (16) 
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where  higher  order  terms  in  a  have  been  neglected  and  <<:>>  represents  the 
second  time  average  over  the  v  ^  time  scale.  This  change  in  the  particle 
energy  is  simply  the  increase  or  decrease  of  the  sloshing  energy  of  the 
electrons  in  the  wave  as  the  wave  amplitude  changes  and  is  entirely  re¬ 
versible. 

The  previous  discussion  of  the  small  Larmor  radius  particles,  which  are 
localized  near  the  0-points  can  be  easily  generalized  to  arbitrary  Larmor 
radius.  Since  the  change  in  the  constant  H  curves  in  Fig.  la  is  much 
slower  than  the  rapid  cyclotron  motion  around  the  curves,  the  area  within  a 
particle  orbit  must  be  preserved,  i.e. 


j  =  1  e  de  (17) 

i 

is  constant.  In  the  previous  calculation  for  electrons  near  the  0- point,  the 
area  within  the  orbit  could  only  be  preserved  if  the  electron  followed  the 
0-point.  To  demonstrate  that  no  irreversible  electron  heating  can  occur  for 
ct  <  1  and  arbitrary  Larmor  radius,  we  allow  the  wave  amplitude  a  to  increase 
to  some  maximum  amplitude  and  then  decrease  to  zero.  Since  the  area  within 
the  particle  orbit  is  the  same  in  the  initial  and  final  states  and  the  energy 
is  a  single  valued  function  of  the  area,  no  net  change  in  the  particle  energy 
can  occur  so  all  electron  "heating"  is  reversible. 

B.  Large  Amplitude  Wave;  a  >>  1 

When  v  ^  0,  the  islands  shown  in  Fig.  lb  for  a  >  1  move  toward  positive 
9  with  an  approximate  velocity  9  %  v,  the  phase  velocity  of  the  wave.  More 
specifically,  an  island  is  formed  at  0  =  -a  with  zero  amplitude.  The  x-point 
propagates  with  a  velocity 
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0X=v(l+a)  (J8) 

s 

and  the  O-point  with  a  velocity 

0°  -  v(l  -  a"1)  .  09) 

s 

Thus,  the  island  half  width  Aft  increases  as  the  island  propagates  at  a  rate 

A0  £  2v/ce  .  (20) 

The  island  reaches  maximum  amplitude  at  G  =  0  and  then  shrinks  at  the  same 

rate  and  disappears  at  0=ot.  The  lifetime  x  of  a  given  island  is  therefore 

s 

r  2-  2ct/v  ,  (21) 

s 

The  motion  of  electrons  which  are  far  from  the  separatrices  in  Fig.  lb 
can  be  calculated  in  a  fashion  analagous  to  the  weak  field  limit  a  <  1.  For 
example,  an  electron  within  an  island  in  Fig.  lb  must  move  toward  positive  0 
to  preserve  the  J  invariant  of  its  orbit.  However,  since  the  island  shrinks 
and  eventually  disappears  as  it  propagates  toward  0=a  ,  all  electrons  within 
a  given  island  must  eventually  cross  the  separatrix  of  the  island.  As  the 
electrons  intercept  the  separatrix,  the  J  invariant  is  broken  sinre  the  period 
of  the  particle  orbit  on  the  separatrix  is  infinite  and  therefore  nu  longer  small 

compared  with  v  The  detrapping  or  conversely  the  trapping  of  electrons  is 

2 

therefore  an  irreversible  process.  All  electrons  with  energy  H  <  a  /2 
eventually  participate  in  this  trapping-detrapping  process  and  we  would  there¬ 
fore  expect  the  electrons  to  eventually  populate  this  entire  phase  space.  De¬ 
pending  on  the  initial  electron  temperature,  this  strong  wave-particle  inter¬ 
action  produces  substantial  bulk  electron  heating.  In  unnormalized  units,  the 
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effective  electron  temperature  becomes 

T  £  m  cV/B2  ,  (22) 

e  y 

i.e.,  the  electron  thermal  velocity  is  simply  cE^/B. 

Electron  trapping  by  a  large  amplitude  perpendicularly  propagating  wave 

has  been  observed  previously  in  computer  simulations  of  the  electron  beam 

12 

cyclotron  instability.  The  qualitative  features  of  the  trapping-detrapping 

process  for  a  >>  1  were  observed.  In  these  simulations,  however,  the  wave 
frequency  and  electron  cyclotron  frequency  were  comparable  (v  'v  1)  and  J  was 
not  a  constant  of  the  motion.  We  would  expect,  therefore,  that  electron  trapping 
would  be  much  less  complete  than  in  the  case  v  <<  1. 


III.  MULTIPLE  WAVES 


The  invariance  of  J  for  the  motion  of  electrons  in  a  low  frequency  wave 
w  «fte  was  essential  in  deriving  the  stochasticity  threshold  in  Eq.  (1).  As 
long  as  the  low  frequency  assumption  is  satisfied,  the  previous  calculation 
can  be  extended  to  include  an  arbitrary  number  of  waves.  The  equation  of 
motion  for  electrons  in  an  arbitrary  spectrum  of  one-dimensional,  low-frequency 
waves  is  simply 


y  +  y  =  -c  Ey(y,t)/Bfi2 

=  -(c/Bfi2)  l  Ek  cos (ky  -  vfct  +  <J)R) 
k 


(23) 


with  an  energy 


H  =  y2/2  +  y2/2  -  (c/Bfi2)  £  E  sin(ky  -  v  t  +  <)>  )/k 


e  £  k 


(24) 


The  electron  motion  can  again  only  become  stochastic  once  multiple 
stationary  points  are  formed,  which  requires 


3E  / 9y  >  i|  B/c 
y  e 


(25) 


When  this  threshold  is  exceeded  locally,  electron  heating  can  occur.  In  most 
physical  systems  of  interest,  precise  wave  phase  information  is  not  available 
and  a  threshold  based  on  averaged  stochastical  properties  of  the  wave  turbu¬ 
lence  is  desireable.  If  the  phases  <}>k  of  the  individual  waves  are  random, 
we  simply  square  Eq.  (25)  and  average  over  the  random  phases  to  obtain  the 


stochasticity  condition 


14 


y  I cE  /B I  2  k 2/ttl  >  1  .  (46 ) 

,  k  e 

k 

Particle  motion  in  random-phased,  perpendicularly  propagating  waves  has 
been  discussed  by  a  number  of  authors  in  connection  with  turbulence  theories 
based  on  the  renormalization  of  particle  propagators.  Dum  and  Dupree  investi¬ 
gated  particle  motion  in  a  general  spectrum  of  low  frequency  fluctuations.^ 

In  their  formalism,  electrons  can  undergo  spatial  diffusion  in  a  one¬ 
dimensional  spectrum  of  waves  once  the  following  threshold  is  exceeded, 

l  |  cE  /B  j 2  kV/w2  >  1  ,  (27) 

k  k 

where  p  =  v  /fl  is  the  electron  Larmor  radius.  In  contrast,  the  present 
e  e 

calculation  demonstrates  that  the  spatial  excursion  of  electrons  is  always 

bounded  since  the  constant  H  curves  are  always  closed.  Spatial  and  velocity 

diffusion  are  only  possible  in  the  limited  sense  that  the  electron  position 

and  velocity  can  scatter  over  a  bounded  region  of  phase  space  [for  example, 

2 

H<  a  in  Fig.  lb]  and  only  when  the  threshold  given  in  Eq .  (26)  is  exceeded. 

A  more  detailed  discussion  of  the  absence  of  spatial  diffusion  in  this  one¬ 
dimensional  model  is  presented  in  Appendix  A. 
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IV.  APPLICATION  TO  THE  LOWER  HYBRID  DRIFT  INSTABILITY 


The  basic  properties  of  the  lower-hybrid-drift  instability  have  been 

briefly  mentioned  in  the  Introduction  and  have  been  extensively  discussed  in 
1-3 

the  literature.  The  instability  essentially  causes  the  plasma  density 

profile  to  become  fluted  perpendicular  to  B.  Below  the  thresholds  given 
in  F.q.  (2)  in  the  case  of  a  single  wave  and  Eq .  (24)  for  the  case  of  manv 
waves,  the  electrons  simply  undergo  coherent  oscillations  in  the  wave  and  the 
lower-hybrid-drift  wave  cannot  irreversibly  exchange  energy  or  momentum  with 
the  electrons.  Thus,  if  the  instability  saturates  before  this  threshold  is 
exceeded,  the  lower-hybrid-drift  instability  cannot  cause  anomalous  diffusion 
of  an  inhomogeneous  plasma.  In  this  limit,  the  plasma  essentuallv  evolves 
to  a  complicated  fluted  state  but  the  mechanism  by  which  this  fluted  state 
can  evolve  to  a  diffuse  profile  must  involve  new  physics  not  yet  incorporated 
into  the  present  model. 

The  distinction  between  the  evolution  of  an  inhomogeneous  plasma  to  a 
fluted  state  and  the  evolution  to  a  more  diffuse  profile  has  not  been 
adequately  addressed  within  the  literature,  possibly  because  it  has  been 
assumed  that  once  the  flutes  develop,  theevolurion  to  a  broader 
smooth  profile  was  an  inevitable  process.  In  computer  simulations  of  the 
lower-hybrid-drift  instability,  for  example,  the  formation  of  these  flutes  has 
been  observed.  The  density  profile  is  then  calculated  by  averaging  over 
these  flutes.  This  averaging  process  can  produce  an  apparent  broadening  of 
the  plasma  profile  even  though  no  diffusion  has  actually  occurred.  To 
illustrate  this  point,  we  consider  a  periodic  displacement  of  the  density 
profile 


n(x,y)  =  nQ  [ x  +  Ax(y) ]  . 


i 


where  x(v)  Is  the  displacement.  ('learlv,  iu>  diffusion  has  taken  place  since 
for  a  given  v,  the  plasma  profile  is  unchanged.  Yet  when  n(x.v)  is  averaged 
over  v,  we  find 

2  2  2 

•n-  <n(x,y)->  =  n  (x)  +  <Ax‘  •  fi  ^  n  I  .\nC  .  (29) 

v  o  v  o 


Tile  profile  has  effectively  diffused  as  shown  in  Fig.  3.  This  fake  "diffusion1 

is  analogous  to  tiie  electron  "heating"  associated  witli  the  coherent  sloshing 

of  electrons  in  the  wave.  In  computer  simulations  of  this  instability,  very 

little  real  diffusion  is  observed  in  t lie  weak  drift  regime  v  ,v.  while  in 

d  i  t 

the  stronger  drift  regime  v . .  >> v.  real  diffusion  is  observed.2  We  now  con- 

d  i  1 

sider  under  what  conditions  the  thresholds  in  Eqs.  ( 1’ )  and  (26)  can  he 

exceeded,  allowing  transport  to  take  place. 

Three  basic  saturation  mechanisms  for  the  lower-hvbr id-dr  if t  instability 

have  been  discussed:  (1)  trapping  or  flattening  of  the  ion  velocity  distrihu- 
6  13 

tion  function;  ’  (2)  depletion  of  the  free  energy  available  to  drive  the 

4  6  13 

instability;  ’  *  and  (3)  electron  resonance  broadening.  The  1 ower-hvbr id- 
drift  instability  is  basically  a  negative  energy  wave  driven  bv  ion  l.andau 
damping.  Flattening  of  the  ion  distribution  either  hv  quasil  inear  relaxation 
or  trapping  can  therefore  lead  to  saturation  bv  quenching  the  ion  l.andau 
damping. 

The  energy  source  of  the  lower-hybrid-drift  mode  is  the  drift  energy  pro- 

2 

duced  by  the  local  gradients,  nmv^/2.  It  was  previously  predicted  that  the 
wave  energy 


W  =  w(3t/3<1.)(E2/8n) 
h  v 


(30) 


=  2(1  +co2  /Q2)(E2/8ti) 
pe  e  y 

4 

could  not  exceed  the  particle  drift  energy  or  l ha l  the  electric  field 


fluctuation  must  satisfy 
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(cE  /B)2  -  v  2  /[2(1  +  2/J  )]  .  (31) 

y  di  e  pe 

In  a  finite  beta  system,  of  course,  the  drift  energy  and  magnetic  energy  are 
linked  so  that  the  magnetic  energy  also  changes  as  the  local  gradients  change}5 
In  a  straight  e-pinch,  the  change  in  the  magnetic  energy  AW  as  the  gradient 
relaxes  is  approx ima t e I v  (see  Appendix  B) 

AWB  =  n<Te+Ti)  R(B)  (32a) 

wi  tli 

g(p)=i(i+fi)1/2-n/i(i+B)1/2+n.  (32b) 

> 

where  f  =  RTrn(T  +T.)B“  is  a  representative  value  of  the  local  plasma  beta, 
e  i  ' 

1/9 

This  magnetic  free  energy  exceeds  the  drift  ene.rgv  as  long  as  v/v  <  (m./m  ) 

d  i  i  i  e 

so  the  wave  energy  hound  becomes 

(<VR)2  '  VL  R(P)/n+f2/ui2e]  ,  (33) 

> 

where  v“  =  (T  +T.)/m  .  In  a  system  where  the  particle  drifts  support  a  reversed 
es  e  i  e 

magnetic  field,  such  as  a  reversed  field  9-pinch,  the  gradients  cannot  be 
relaxed  without  dissipating  all  the  magnetic  energy  in  the  system,  the  free 
energy  available  to  drive  the  Instability  is  effectively  infinite  and  there  is 
no  wave  energy  bound}R 

Resonance  broadening  lias  also  been  proposed  as  a  saturation  mechanism  tor 

14 

the  lower-hybrid-drift  instability.  Once  the  fluctuating  electric  fields 
exceed  a  critical  amplitude,  the  electrons,  which  are  initially  nonresonant, 
can  effectively  resonantly  interact  with  the  waves.  The  resulting  "Landau" 
damping  of  the  wave  energy  stabilizes  the  wave  spectrum.  The  threshold  at 
which  this  transition  takes  place  is,  of  course,  the  subject  of  this  paper 
and  is  given  by  either  Eq .  (2)  or  Eq .  (2b)  in  the  case  of  single  or  multiple 
waves,  respectively.  These  thresholds  are  substantially  higher  than  those 
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predicted  on  the  basis  of  resonance  broadening  theory  | Eq .  (27)). 

Above  the  stochasticlty  threshold,  strong  electron  heating  takes  place,  thus 

dissipating  the  wave  energy  and  stabilizing  the  lower  hybrid  drift  wave.  If 

ion  trapping  or  depletion  of  the  free  energy  do  not  occur  first,  the  lower- 

hybrid-drift  instability  will  saturate  at  an  amplitude  given  by  either  Eq.  (2) 

or  Eq.  (26).  {n  the  previous  application  of  the  resonance  broadening  theory 

to  the  lower-hybrid-drift  instability,  substantially  lower  thresholds  were 

14 

predicted  when  electron  VB  drifts  were  included.  The  question  of  the 
validity  of  the  resonance  broadening  theory  in  the  presence  of  the  VB  drifts 
cannot  be  addressed  on  the  basis  of  the  present  calculation.  The  reduction 
of  the  stochasticity  threshold  by  the  VB  drift,  however,  seems  physically 
reasonable,  and  the  inclusion  of  these  effects  would  be  an  important  extension 
of  the  present  work. 

We  now  consider  under  what  conditions  saturation  of  the  lower-hvbrid- 

drift  instability  can  occur  by  stochastic  electron  heating.  The  growth  rate 

of  the  lower-hybrid-drift  instability  peaks  in  the  range  kpgs  -  1-2,  where 
2  2  2 

p  -v  /fl  ,  so  the  electron  stochastiritv  threshold  in  Eq .  (2)  becomes 
es  es  e  - 

e*/T.  %  (kp  )~2  =  -25  -  1  .  (34) 

l  es 

In  the  weak  drift  limit  v^/v^  <<  1>  the  phase  velocity  of  the  wave  is  Vp^  'v 
vdi<<  V  The  saturation  of  the  instability  can  occur  by  flattening  the  ion 
velocity  distribution  or  ion  trapping  at  an  amplitude 

e*  =  mivph/2  =  Vdi72  C<  Ti  (35) 

13 

even  in  the  absence  of  the  wave  energy  bound.  This  amplitude  is  below  the 
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threshold  for  electron  stochasticity  given  in 


is  possible  in  the  weak  drift  regime.  In  the 

trapping  occurs  when  e<t>  -c  m.v“,  /I  >  m.v^/2  or^ 

l  ph  11 


Eq  .  (  14  )  so 
strong  drift 


no  electron 
regime  vdj 


<i i  £  f  usion 
^  v  ^  ,  ion 


e<f>/T^  >  1 


06) 


so  the  electron  stochasticity  and  ion  trapping  conditions  are  comparable. 

In  the  case  of  a  reversed  magnetic  field  in  the  strong  drift  regime  both 
ion  trapping  and  electron  stochasticity  should  cause  saturation  of  the  in¬ 
stability  and  stromt  electron  heating  and  dit  fusion  »  an  occur.  For  k  t  1, 

the  wave  energy  bound  in  Eq.  (33)  for  the  straight  pinch  can  he  rewritten  as 

e<J>/T .  <  g(f  )  .  (37) 

In  a  high  8  configuration  this  hound  is  comparable  to  the  ion  tripping  and  electron 
stochasticity  conditions  so  that  diffusion  and  heating  will  occur.  Tn  a  low 
g  configuration  the  energy  bound  is  smaller  than  either  of  these  conditions 
and  no  diffusion  is  possible. 

Two-dimensional  particle  simulations  of  the  lower-hybrid-drift 
instability  have  been  carried  out  in  both  straight  and  reversed  magnetic 
fields. These  runs  were  carried  out  for  relatively  strong  drift  1  10 

with  relatively  large  8 'vO. 25-1.0  and  artificial  mass  ratios.  As  a  consequence, 
the  wave  energy  bounds  in  Eqs.  (31)  and  (37),  the  ion  trapping  threshold  in 
Eq.  (36)  and  the  electron  stochasticity  threshold  In  Eq .  (34)  are  all 
comparable  even  for  the  straight  pinch.  Electron  heating  should  therefore  be 
expected,  especially  in  the  strong  drift  regime  v^/v^  >>  1 •  In  the  case  of 
v^/vj^l  fluting  of  the  density  profile  is  observed  in  the  simulations  hut  the 
mode  saturates  with  very  little  actual  broadening  of  the  profile. ^  In  the 
strong  drift  regime  both  strong  electron  heating  and  diffusion  are  observed. 


as  expected. 
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V .  SUMMARY  AND  CONCLUSIONS 

Electron  motion  in  a  low-frequency  (ui  <<  s;^)  wave(s)  propagating  per¬ 
pendicular  to  a  uniform  magnetic  field  is  studied.  Both  monochromatic  and 
broadband  electric  field  spectra  have  been  considered  although  the  investigation 
is  strictly  limited  to  one-dimensional  spectra.  Below  a  threshold  amplitude 
given  by 


<i  =  k  c  E  /Bil  ~  1  (37) 

y  y  e  '  ' 

for  the  case  of  a  monochromatic  wave  and  the  corresponding  generalization 

l  Ik  cE  / B:j  | 2  -  1  (38) 

k  K  e 

for  a  broadband  spectrum  of  waves,  electron  heating  and  momentum  transfer 
are  strictly  reversible.  Above  these  thresholds,  the  Larmor  motion  of  an 
electron  in  a  uniform  magnetic  field  is  strongly  modified  and  trapping  of  the 
electrons  by  the  wave  can  occur.  The  electrons  remain  trapped  for  a  time  t 
given  by 

r  '  k  cL  /Bui;:  =  u/w  a  .  (39) 

tr  y  y  e  6 

The  process  of  trapping  and  detrapping  allows  an  irreversible  exchange  of  energy 
between  the  electrons  and  the  wave  and,  in  particular,  causes  bulk  "heating" 
of  electrons  to  a  mean  velocity 

v  v  v  =  cE  /B  . 
l  h  y 


(AO) 


This  crapping-detrapping  phenomena  has  been  previously  observed  in  computer 

simulations  of  the  electron  beam  cyclotron  instability,  although  the  trapping 

12 

is  much  less  complete  in  this  limit  because  oj  'v  Si  . 

e 

The  implications  of  these  results  for  the  nonlinear  evolution  of  the  lower- 
hybrid-drift  instability  and  associated  transport  have  been  explored.  The  growth 
rate  of  the  lower-hybrid-drift  instability  peaks  for  kp  =  k(T./m  )1/,2/ii  -  1-2 

GS  1  G  G 

so  the  threshold  in  Eq.  (37)  becomes 

e4>/T.  =  .25  -  1  .  (41) 

Below  this  threshold  no  irreversible  electron  heating  or  momentum  transfer  is 

possible  and  the  lower-hybrid-drift  instability  cannot  cause  diffusion  of 

electrons  in  an  inhomogeneous  plasma.  The  apparent  diffusion  which  is  found 

from  quasilinear  theory  in  this  limit  corresponds  to  the  formation  of  flutes 

or  ripples  in  the  plasma  profile — not  to  the  evolution  to  a  more  diffuse 

profile.  In  the  weak  drift  regime  vcjj_<vi>  the  lower-hybrid-drift  instability 

saturates  below  the  threshold  given  in  Eq.  (41)  by  trapping  ions  or  depleting 

13 

the  source  of  free  energy,  and  transport  by  the  lower  hybrid  drift 
instability  is  not  possible  (at  least  within  the  limitations  of  the  present 
model) . 

In  the  strong  drift  regime  v^>v^,  the  lower  hybrid  drift  instability 
grows  until  the  threshold  in  Eq.  (41)  is  exceeded.  The  rapid  onset  of  strong 
electron  heating  above  this  threshold  prevents  further  amplifications  of  the 
wave  and  thus  saturation  of  the  instability  occurs  for  e<f>/T^<l.  Electron 
heating,  transport,  and  magnetic  energy  dissipation  accompany  the  saturation 
of  the  mode  in  this  limit ^  and  have  been  observed  in  recent  computer  simula¬ 
tions  of  the  lower-hybrod-dr if t  instability  in  a  reversed  magnetic  field. ^ 
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It  should  be  emphasized  that  these  conclusions  are  strictly  based  on  the 
present  limited  model:  (1)  a  one-dimensional  wave  spectrum;  (2)  k((  =  k*jB/|B|  =  0; 
and  (3)  a  uniform  magnetic  field.  The  relaxation  of  these  assumptions  should 
lead  to  a  reduction  of  the  threshold  given  in  Eq .  (31)  and  allow  the  lower- 
hybrid-drift  instability  to  cause  anomalous  transport  in  the  weak  drift  regimes. 
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APPENDIX  A 


We  now  consider  in  more  detail,  the  electron  dynamics  in  a  one-dimensional 
spectrum  of  low  frequent'-  (10  <<  Wg)  waves.  The  electron  equation  of  motion  is 
given  in  Eq.  (23), 


v  +  y  =  -cEy(y,  t)/BS2e  I  -e(y,t) 


(Al) 


We  have  shown  that  the  electron  motion  is  bounded  and  y  therefore  remains 
finite  for  all  time.  For  simplicity,  we  limit  our  investigation  to  the  case 
where  the  stochasticity  threshold  in  Eq.  (26)  is  not  exceeded  and  the  electron 
Larmor  radius  is  small  compared  to  the  wavelength  of  the  fluctuations.  In  this 
case,  e(y,t)  can  be  expanded  around  y  =  0, 


y  +  y (1  +  3e/3y)  =  -e(0,t) 


(A2) 


The  electric  field  fluctuations  cause  a  frequency  shift  of  the  electron  gyro- 


motion  5W(t)  =  (3e(y,t)/3y) 


y=0 


The  solution  of  Eq.  (A2)  is 


y(t)  -  -e (t)  +  pe  cos(t  +  <p  +  6<j>), 


(A3) 


where  is  the  electron  Larmor  radius,  4>  is  the  initial  phase  angle  of  the 
electron  in  its  gyromotion  and 


<5<j)  (t)  = 


ft 

di 

■’o 


6W(t) 


(A4) 


is  the  phase  shift  induced  by  the  fluctuating  electric  field.  Averaging  y(t) 
over  an  ensemble  of  fluctuations, 
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<  y (t)>  =  p  <-  cos(t  +  +  <H‘ ) 


P  coi'.(t  +  <))  exp l -  4'(  t )  / 2 ]  , 


*(t)  ■=  j  dT1d,2<6j(-rJ)6S2(:2):  , 


<>  derotes  an  average  over  the  ensemble  and  the  cumulanL  expansion  has  been 
used  to  evaluate  <  exp(i54>)  >  -  txp(-4/2).22  T)ie  electron  diffusion  can  similarly 
be  calculated. 


<  6y2>  =  <  (y  -  <  y  ■)  **>  =  t-2>  +  (r  V2)  [  1  -  exp(-  1)  ]  . 


If  the  electric  field  fluctuations  have  a  finite  correlation  time  t  ,  then  for 

c 

t  >  Tc, 


<6y2>  =  <  t:“>  +  (p72)[l  -  exp(-4)J 


where 


i  =  (2ir)  1  dk  do)  5  (to)  I  Si;,  I2  . 

k,u 1 


The  first  term  in  Eq .  (AS)  corresponds  to  the  mean  excursion  of  the  electron 
associated  with  its  polarization  drift.  The  second  term  arised  because  ol.  the 
frequency  shift  of  the  electron  Larmor  motion.  The  phase  of  the  electrons 
along  their  orbits  becomes  uncertain  and  is  reflected  in  an  uncertainty  in  their 
position.  Over  short  times  this  appears  as  a  diffusion 


2  2  • 

<  6y  >  %  p  $  t/2 


(A10) 


of  the  particle  position  [compare  with  Eq .  (7)  in  Kef.  14  in  the  long  wavelength 


I 


limit].  However,  over  longer  times  it  >  1,  the  maximum  uncertainty  in  the  particle 
2 

position  is  Pe/2,  i.e.,  the  phase  angle  of  the  electron  in  its  Larmor  orbit  is 

completely  unknown.  Thus,  spatial  diffusion  is  possible  below  the  threshold  given 

in  Eq  .  (26)  for  i  finite  hut  is  limited  to  the  electron  Larmor  radius, 

c 

This  result  is  consistent  with  the  previous  discussion  of  the  preservation  of  the 
J  invariant  in  Fig.  la.  The  phase  of  the  electron  along  the  curves  shown  in 
Fig.  1  can  become  uncertain  but  the  area  within  the  curve  is  still  preserved  and 
thus  diffusion  is  strictly  limited. 

It  should  also  be  noted  that  even  the  limited  diffusion  described  by  Eq .  (A8) 
is  produced  only  by  the  zero  frequency  component  of  the  wave  spectrum  [see  Eq .  (A9)]. 
The  frequency  shift  of  the  to  ^  0  waves  is  periodic  and  causes  no  net  phase  shift 
of  the  electrons  in  their  Larmor  motion.  In  a  spectrum  of  high  frequency  waves 
no  diffusion  occurs  unless  the  threshold  in  Eq.  (26)  is  exceeded.  The  "bootstrap" 
diffusion  calculated  by  Dum-Dupree  when  the  threshold  in  Eq.  (27)  is  exceeded  does 
not  take  place  because  the  excursion  of  electrons  in  the  fluctuating  electric 
fields  is  strictly  bounded  (there  is  no  orbit  secularity) . 
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APPENDIX  B 

The  lower  hybrid  drift  instability  is  driven  hv  the  drift  energy  of  the 

2 

electrons,  n  in  v~/2.  However,  in  a  finite  8  plasma  the  magnetic  field  and 
e  d 

the  plasma  currents  are  linked  so  the  plasma  currents  can  not  be  relaxed 

without  changing  the  magnetic  energy.  When  the  plasma  currents  support  a 

reversed  magnetic  field,  the  plasma  current  can  not  be  dissipated  without  dis-  i 

sipating  all  the  magnetic  energy  in  the  system  so  there  is  effectively  no 

1 6 

energy  bound  for  the  lower  hybrid  drift  instability.  We  now  consider  the 
change  in  magnetic  energy  as  the  plasma  profile  broadens  in  the  straight  0 
pinch. 

For  simplicity,  the  plasma  profile  is  modeled  by  the  simple  step  shown 
in  Fig.  4  (the  profile  is  taken  to  he  symmetric  around  x=0) .  An  initial 
plasma  of  density  n.,  temperature  and  in  a  magnetic  field  is  supported 
by  an  external  magnetic  field  B^.  The  initial  plasma  width  L*  is  allowed  to 
increase.  Conducting  boundaries  at  X=Lo  prevent  flux  from  entering  or  leaving 
the  system  so  Lite  total  flux 

A  =  B^L*  4  B  (L  -I  ‘)  (Bl) 

it  o  o  i 

is  constant.  The  total  number  density  N = n1L^  is  also  conserved.  The  change 
in  the  particle  drift  energy  is  ignored  since  this  can  only  increase  the  free  I 
energy  available  to  drive  the  lower  hybrid  drift  instability. 

Local  pressure  balance 

=  B^/8n  +  nT  (B2) 


will  also  be  maintained  during  the  evolution  of  the  profile.  Flux  and  particle 
conservation  are  invoked  to  eliminate  B^  and  n  from  (B2)  to  obtain  an  expression 
for  T, 


T/T  = 


[2BoLrLi(VBi)  I/rLi(Bo +Bi)] " 1  • 


(B3) 
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When  L,  >>  L*  T/T1  =  IB  /(B  +b!)  or  the  increase  in  temperature  AT  =  T-T1  is 
11  oo! 

given  by 

AT  =  T-Tl  =  T1|(l+B1)1/2-1  l/Kl+eV^+l]  (BA) 

i  i  i  i2 

where  6  = 8im  T  /B^  is  the  initial  plasma  beta.  In  a  low  6  system  the 
temperature  increase  is  given  by 

AT1  i  Tl  BX/2  <<  T1  .  (B5) 

In  a  high  B  system,  however,  the  plasma  heating  can  be  substantia]  with 

AT1  ^  Tt . 

This  increase  in  the  plasma  thermal  energy  is  matched  by  a  corresponding 
reduction  in  the  magnetic  energy,  i.e.,  as  the  sheath  broadens,  the  system 
evolves  to  a  lower  magnetic  energy  state.  This  energy  must  be  included  in 
the  free  energy  available  to  drive  the  lower  hybrid  drift  mode  since  the 
plasma  currents  can  not  be  relaxed  until  all  of  the  magnetic  energy  has  been 
dissipated.  Although  the  expression  for  AT  in  Eq.  (BA)  was  derived  for  the 
rather  simplified  step  model  in  Fig.  A,  the  qualitative  scaling  with  B 

is  valid  for  a  more  general  profile  and  will  be  discussed  in  a  subsequent 

...  .  15 

publication. 
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FIGURE  CAPTIONS 


Constant  H  curves  are  shown  in  the  8-0  phase  plane  for  a  =  .5, 

$  =  0  in  (a)  and  «  =  8tt  ,  <ji  =  it / 2  in  (h)  . 

The  positions  of  the  stationary  points  0g  =  a  cos  (8^  +  <}>)  are  shown 
graphically  for  the  parameters  of  Fig.  lb. 

The  pseudo-diffusion  <n  [x+Ax(v)]>  =n  (x)  + (l/2)n"<Ax  >  which 

O  y  o  o  y 

results  from  the  periodic  displacement  Ax(y)  of  an  initial 

plasma  profile  n  (x)  is  shown.  The  bracket  <>  denotes  a  spatial 

o  y 

average  over  y. 

The  8-pinch  profile  is  modeled  by  a  simple  step  in  the  magnetic 
field  profile.  The  profile  is  symmetric  around  x=0. 


